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Abstract 

We study the scalar and tensor integrals associated with the pentabox topology: 
the class of two-loop box integrals with seven propagators - five in one loop and three 
in the other. We focus on the case where the external legs are light-like and use 
integration-by-parts identities to express the scalar integral in terms of two master- 
topology integrals and present an explicit analytic expression for the pentabox scalar 
integral as a series expansion in e = (4 — D)/2. We also give an algorithm based on 
integration by parts for relating the generic tensor integrals to the same two master 
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1 Introduction 



Scattering processes are one of the most important sources of information on short dis- 
tance physics and have played a vital role in establishing the properties of the fundamental 
interactions of nature. At the one-loop level, box integrals are a key ingredient in the 
comparison of perturbative predictions with experimental data, such as wide-angle jet pro- 
duction in hadron collisions or Bhabha scattering in electron-positron annihilation. Recent 
improvements of experimental measurements demand even more precise theoretical predic- 
tions and there is significant interest in determining 2 — > 2 scattering rates at the two-loop 
order. Achieving this goal requires the evaluation of two-loop graphs such as the planar 
double-box graph JI|, || , the non-planar double-box graph or some one- loop box integrals 
with bubble insertions on one of the propagators 0]. 

In addition, we also need to study the particular type of planar box graph shown in 
Fig. [I]. Owing to the particular shape of this diagram, we name this topology pentabox. 
It is the purpose of this paper to provide analytic expressions for the scalar pentabox with 
unit powers of the propagators using dimensional regularisation for the light-like (or on- 
shell) case, y\ = y\ = v\ = v\ = 0, in terms of generalised polylogarithms. This type of 
integral is much simpler than the planar double-box and the non-planar double-box []3| 
because of the presence of the triangle sub-graph which allows the use of integration-by- 
parts identities || to relate the most general scalar integral to two types of simpler integrals 
with fewer propagators, which we denote master topologies. We give expressions for these 
simpler integrals and use them to derive an analytic expression for the pentabox graph 
with unit propagators as a series expansion in e = (4 — D)/2. We also give algorithms for 
evaluating all the associated tensor integrals. These are directly related to scalar integrals 
with higher powers of the propagators in higher dimension. As with the scalar integral, the 
integration-by-parts identities allow the tensor integrals to be obtained in terms of the two 
master-topology integrals. 

Our paper is organised as follows. In Section |^ we establish the basic notation we will use 
for two-loop integrals. Analytic expressions for the scalar pentabox integral, together with 
algorithms for reducing the integral to the simpler master topologies, are given in Section |3], 
while the tensor integrals are described in Section f|. Details of the master topologies and 
algorithms relating them are given in Section |5|. Explicit expansions in e = (4 — D)/2 
are computed in terms of generalised polylogarithm functions. Useful relations amongst 
the polylogarithms are collected in Appendix Finally, our findings are summarized in 
Section ||. 
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Figure 1: The pentabox topology. The propagators are labelled according to Eq. ($.9J 
and are each raised to the V{ power. 



2 Notation 



We denote the generic two-loop integral in D dimensions with n propagators Ai raised 
to arbitrary powers as 

d &\ f d [lj k-^ - k<2 •) k\ k\ ) ^2 ' * * *] 



VK 



D/2 



in 



D/2 



Ai ■ ■ ■ A£ 



J D ({^j}; {Qi}) [l; ^i j k 2 j • • •] 

(2.1) 

where the external momentum scales are indicated by {Q 2 }. When the numerator in Eq. Q2.1p 
is unity, J D ({^i}; {Ql}) [1] = J D ({^i}; {Qf}), we have the scalar integral that forms the ba- 
sic underlying skeleton of the graph. The tensor integrals J D {{^i}] {Qf}) [k\ ; . . .], which 
naturally arise in Feynman diagrams, are usually related to combinations of the scalar in- 
tegral with different powers of propagators and/or different values of D [0. Evaluation of 
two-loop matrix elements for physical processes requires a knowledge of both tensor and 
scalar integrals. 

In this paper we focus on the specific two-loop box graph shown in Fig. [I], whose analytic 
expression is 



J (f i, v 2 , v 3 , 1/4, "5, "a, v i\ s, t) 



d D h f d D k 2 



1 



in o/2j m D/2 A^A^AfA^AfAfAy' 
together with the associated tensor integrals, where the propagators are given by 

,2 



(2.2) 



A 1 = kf + zO, 
A 2 = {k 1 +p 1 ) 2 + i0, 
A3 = (k 1 +p 1 +p 2 ) 2 + 2O, 
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A 4 = (h + Pl + p 2 + p 3 ) 2 + iO, (2.3) 

A = (A; 2 +Pi+P2+P3) 2 + «0, 

A 6 = k 2 2 + iO, 

A 7 = (h-hy + iO. 

With this momentum assignment we see that the loop momentum k\ circulates in five of the 
propagators while k 2 circulates in only three. All of the external momenta are in-going and 
are taken to be light-like: p\ = p\ = p\ = pi = 0. The only momentum scales present in 
the problem are therefore the usual Mandelstam variables s = (p\ + p 2 ) 2 arid t = (p2 + P3) 2 , 
together with u = —s — t. 



3 The scalar pentabox integral 

The generic scalar pentabox integral, J D (ui, z/ 2 , u 3 , i/ 4 , v 5 , z/ 6 , u 7 ; s, t), with arbitrary pow- 
ers of propagators in D dimension is always reducible to a finite set of simpler master- 
topology integrals. In fact, applying the integration- by-parts method || to the following 
two integrals 

rdPh rd D k 2 d i(k 2 - P ,r- 

J m D ' 2 J m°l 2 dk^ A^A^A^A^A^A^A^ ' 1 } 

we find that 

(D - 2z/ 5 - u 6 - vi) J D = (z/ 6 6 + 5 - + u 7 7+5~ - v 7 7 + 4~) J D , (3.2) 

(D - u 5 - 2z/ 6 - v 7 ) J D = (y b h+Q- + v 7 7 + 6~ -v 7 7 + l~) J D , (3.3) 

where J D = J D (ux, v 2y ^3, ^4, ^5, V7', s, t) and the symbols i + and i _ are shorthand notation 
for raising and lowering of powers of propagator % 

i±J D (...,z/,,...) = J D (...,z/,±l,...). (3.4) 

As usual, each raising operator is always accompanied by a factor of ^ so that it is impossible 
to raise the power of the propagator if it is not already present, i.e. ^ 7^ 0. 

By repeated application of Eq. ( |3.2| ), we can reduce either of or u 5 to zero. Similarly, 
by applying Eq. ( |3.3|) we can lower (and eventually eliminate) the power of either v\ or uq. 



This produces a pinched graph of the form 

J D (^1,^2,^3,^4,0, 0,v 7 ;s,t) = 0, (3.5) 

J D (^1,^2,^3,0, 1/5,0, u 7 ;s,t) = A D (v 1 ,v 2 ,v 3 ,v 5 ,i> 7 ;s,t) , (3.6) 

J D (0, z/ 2 ,z/ 3 ,z/ 4 ,0,z/ 6 ,z/ 7 ;s,t) = B D (u 2 ,u 3 ,u 4 , u 6 ,u 7 ;s,t), (3.7) 

J D (0, z/2,^3,0, u 5 ,u 6 ,u 7 ;s,t) = C D (u 2 ,u 3 ,u 5 ,u 6 ,u 7 ;s,t) , (3.8) 
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where, diagrammatically, 



Pi 



A D (u u u 2: u 3 ,u 5: u 7 ;s,t) = 



P2 



Pi 



B D (u 2 , u 3 , z/ 4 , v 6 , v 7 \ s,t) 



P2 — * 




P4 



P3 



P4 



< P3 



Pi > 



C (v 2 ,v 3 ,is 5 ,v 6 ,is 7 ;s,t) 



P2 




< P4 



P3 



A D , B D and C D constitute our set of master-topology integrals. Due to the particular 
symmetry of these diagrams, they satisfy some additional identities: 

- The basic A D and B D integrals are related by exchanging s and t according to 

B D (u 2 ,u 3 ,u A , v & ,vr,s,t) = A D (is 2 ,is 3 ,v 4 ,v 6 ,is 7 ;t,s), (3.9) 
and A D has the additional symmetries 

A D {y 3 ,v 2 ,v x ,v^v 7 \s,t) = A D (z/i,z/ 2 ,z/ 3 ,z/ 5 ,z/ 7 ;s,t), 

A D (vi, u 2 , v 3 , v 7 , i/ 5 ; s,t) = A D {v x ,v 2 ,v 3 ,v h ,v 7 \s,t). (3.10) 

- Similarly, the basic C D integral has the properties 

C D (v 6 ,v 5 ,is 3 ,is 2 ,is 7 ;t,s) = C D (v 2 ,v 3 ,is 5 ,v 6 ,is 7 ;s,t), 
C D (v 3 ,is 2 ,]s 6 ,is 5 ,is 7 ;t,s) = C D (u 2 ,u 3 ,u 5 ,u 6 ,u 7 ;s,t). (3.11) 



Application: We consider here the specific case of the scalar integral with all the propa- 
gators equal to one. Applying Eqs. ( |3.2|) and ( |3.3| ) we obtain 



j"(l,i,i,l,l,i,W) 



2C D (1, 1, 1, 1, 3; s, t) - 2A D (1, 1, 1, 1, 3; s, t) - A D (1, 1, 1, 2, 2; a, t) 
-2B D (1, 1, 1, 1, 3; s, t) - B D (1, 1, 1, 2, 2; s, t) 



(3.12) 



The explicit form for the master-topology integrals with arbitrary powers of propagators and 
arbitrary dimension are collected in Section ^. Making a series expansion in e = (4 — D)/2, 
we find that the general structure of the scalar integral can be written 



j d (i,i,i,i,i,i,1;M) 



l + 3e) 



r(l -eYT(l + 2e) 



4e 4 (i + e)(i + 2e) 
x —^R(s,t)- 



rfl-3e 



ts 2 



-2c 



-R(t,s) 



(3.13) 



1) In the physical region s > 0, t < we have 

's + t 



R(s,t) 



1 + e + 2C 2 e 2 + 
's + t 



4 Si 



+ 2Cs + 6C2 



+ {l2S 2 , 2 (^)-8S 1 



,3 



s + t 



+ 6C3-51C 4 }e 4 + 4z7re(s,t) | 
3Li 3 



Li 2 



+ 8Si, 2 
s + t 



s + t 



28( 2 Li 2 



s + t 



+ 



s + t 



2 Li, 



s + t 



2 Si 



s + t 



Cs - 2( 2 



+ O (e 5 ) • 
(3.14) 



The generalised polylogarithms (defined in Section |A.1|) that appear in Eq. (|3.14|) , and 
then in Eq. ( p.13] ), have arguments (s + t)/t and (s + t)/s, which are always less than 
unity in this region. The function 0(i, s) assumes values 



e( s ,t) = -e(t, s ) 



(3.15) 



The imaginary parts associated with the prefactor can be simply obtained by recalling 
the +i0 prescription associated with the external kinematic scales 



-s 



-t) 



-2c 



-2c 



s\ e exp (27rie) 

2c 



t 



(3.16) 
(3.17) 
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2) In the other physical region s < 0, t < 0, we use the analytic continuation formulae 
of Section [A.2| , and we find 



R{s,t) = l + e + 2C 2 e 2 + 
+ 3 log 3 



+ |4S 2 , 2 
- 16Li 4 



s + t 
s 



4Li 3 
+ 2tt 2 log 
+ 8S li3 



s + t 
s + t 



4 log 



s + t 



Li-; 



s + t 



s 
s 



+ 4 



log 2 



s + t 
s + t 



+ 4 



2 -lo 



s + t 
s + t 



-4S lj2 
log 



+ 2 log 



s + t 



+ -7T 2 log 2 



S + t 



12[C 3 + 2C 2 ] log 



Li 3 

Li 2 

s + t 



s + t 
s + t 

s 

s 



+ 6Cs + tt' 



+ 1 
+ - lo; 



Si,: 



s + t 



4 /'S + t 



s + t/ 6 



s + t 



4 



log 3 



s + t 



+ 14C 3 + 31C 4 }e 4 + 0(e 5 ). 

(3.18) 



Notice that, in this region, the scalar pentabox integral of Eq. Q3.13| ) has no imaginary 
part, as expected. 



4 The generic tensor integrals 

In this section, we want to deal with the generic tensor integral, J D ({^}; {Q 2 }) \k\ \ • • •], 
using the pentabox just as an example of how to reduce the tensor integral to a set of simpler 
scalar integrals. 

A method to reduce tensor integrals constructing differential operators that change the 
powers of the propagators as well as the dimension of the integral was presented in Ref . pi . 
However, it is in our view simpler to obtain the tensor integrals directly from the Schwinger 
parameterised form of the integral. In fact, in its most general form, the integrand can be 
rewritten as 

^L_ = jv xexv {^, Xi A^ (4.1) 

where 

J fJi T{Vi) Jo 

and 

n 

J2 x i A i = akl + bkl + 2ck 1 - k 2 + 2d - k ± + 2e- k 2 + f, (4.3) 

i=l 

where a, b, c, d^, e M and / are directly readable from the actual graph: a(b) = Y^ x ii where 
the sum runs over the legs in the k\ (k 2 ) loop, and c = J2 x i with the sum running over the 
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legs common to both loops. Specifically, in the pentabox graph we have 



a 
b 
c 
cF 



X\ + X 2 + £3 + £4 + £7 
— £7 

X2P1 + X3P12 + 



X 5Pl23 



f = x 3 s, 



(4.4) 



where pu = Pi + P2 and P123 = Pi + P2 + P3- In all cases, the a, b, c, e M and / are linear 
in the Xj. 



With the change of variables 



kl -> - rill + 



A 9 



where 



and 



cet* - bd» 



V 



cd^ — ae M 
V ' 



V = ab-c 2 



we can diagonalise Eq. ( |4.3|) , so that 



n V O 

Y d x i A i = aKl + -Kl + ± 

i=l a 1 



with 



Q = -ae 2 -bd 2 + 2ce-d + fV. 
The scalar loop integral can be cast in the form 



J D (M; {Ql}) [i] = Jv*j 



d D K\ rd D K, 



exp 



V 



Q 



aKf + -Ki + 

a V 



m D/2 J in D/2 

and the Gaussian integrals over the shifted loop momenta are evaluated to 



J D {M\ {Q-}) [i] =fvxi, 



the integrand X being given by 



2 



T>d/ 



2 ex P 



(4.5) 
(4.6) 

(4.7) 
(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 
(4.13) 



Similarly, the tensor integrals can be easily obtained by using identities such as 
d D K x 



i7T D/2 



K{ exp{aKf) = 0, (4.14) 



d ° Kl K»K\ exp(aK*) = -A gT -L ( 4 - 15 ) 



VK 1 - 



J^m^KW «p(oK?) = ^{sT^ + ^ + grr}-^- (4-16) 
To give a concrete example, we consider the tensor integral associated with fcj* 

Vx X»I. (4.17) 

Recalling the definition (|4.7|) , we see that X^ consists of the ratio of a set of bilinears in Xi 
divided by V. We can therefore absorb the factors of x% into T>x (see Eq. (|4.2| )) by increasing 
the power to which the i-th propagator is raised 

%i -Ui y / ' = (4.18) 



(4.19) 



i> t ) rfa + 1) 

while the factor V can be absorbed into X (see Eq. ( |4.13| )) 

1 1 1 

pD/2 p ==> ' <p(D+2)/2 ' 

so that 1/7 3 acts clS cl dimension increaser 

I =► d +, P =► d-, (4.20) 

where 

d±J D ({i/,}; {Q 2 }) = J D±2 ({i/,}; {<?■}) . (4.21) 



In this way, each Xj in the numerator increases by one the power of the associated propagator, 
and each power of V in the denominator (numerator) increases (decreases) the space-time 
dimension D by two. Schematically we have 

J D (M; {Qi}) [tf] = £ vl J D+2 ({..., ••• , • • •} ; {Q 2 }) [l] , (4.22) 

where the summation runs over the elements of (ce M — M M ) which fix the values of i, j and 
Pk- 
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For the case of the pentabox where a, . . . , f are given in Eq. ( |4.4j ), we have from Eq. ( [4.7[ ) 

x » = -^7^5Pl23 ~ (^5 +X& + X 7 ) (z 2 Pi + X3P12 + S4P123) ^ 33) 

so that Eq. ( [4.17| ) becomes 



-z/ 5 5 + z/ 7 7 + ^ 23 



V 5 5+ + v 6 6+ + z/ 7 7+) + u 3 3+ P 1 2 + z/ 4 4 + ^ 23 )J d+ J D , 

(4.24) 

where we have used the shorthand notation J D = J D ({Vj}; {Qf})- 



By applying the integration-by-parts identities (|3T2|) and ( |3 .3|) we can now write this 
tensor integral directly in terms of the set of simpler integrals. 

For generic four-point integrals, we need tensor integrals with up to four free indices, 
each associated with a Lorentz index of an external leg. Integrals with higher powers of the 
loop momenta are of course possible, but must yield dot products with other momenta when 
the available free Lorentz indices are saturated. In many cases, these dot products can be 
immediately expressed in terms of the propagators and canceled through (see Eq. ( |4.38| )). 

Nevertheless, the procedure previously described can be iterated ad libitum and we can 
express every tensor integrals in terms of scalar ones with increased powers of the propagators 
and dimension D. For example, we have 



J D [kZ] 



J D \KK\ 

J \k\ k 2 ] 
J D [k 2 k 2 ] 

J ^k p k^ k P ^ 

J [kik 1 k 2 ] 

J \k\ k 2 k<2 ] 



/ vx y j, 



2^2] 



Vx 

Vx 
Vx 

Vx 
Vx 

Vx 
Vx 

Vx 



2P y / ' 



x p y v + — <r ) j, 



X »x v X p - — {g^X p + g w X v + g vp X»}j J, 
^X^-^g^ + ^ig^X^ + g^X^ J, 

x p y v y p - ^9 vp x» + ^pW lu y p + g pp y u }) 



y p y v y p - ^pW v y p + g w y v + g up y 1 } ) x, 



(4.25) 
(4.26) 

(4.27) 
(4.28) 

(4.29) 

(4.30) 

(4.31) 
(4.32) 



X^X v X p X a 



47* 



_ {g ^ gP a + gWg ua + gl ,r gVp} 



— {g pv X p X a + g pp X v X a + g pa X v X p + g vp X p, X a 
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+ g va X fl X p + g pa X p X v } j J, (4.33) 

J D [k1klk p x k%\ = fVx I X p X v X p y a - — {g pv X p + g pp X v + g vp X p }y a 
J \ 2"^ 

+ — {g^x v x p + grxw + ^at"*"} 

6c \ 

' {<T<T + W + <T <T} i, (4.34) 



4 p2 



j D \Kk\k p 2 k° 2 \ = Jvx [x*x v yy - ^ j5 g pCT x p x v - ^g^yy 



+ \{g pp x v y a + g vp x p y a + g pa x v y p + g ua x p y p } 

ab c 2 \ 

+ -^sTsT + ^i^sT + f a 9 vp }\ X, (4.35) 



j D [k p k v 2 k p k a 2 ] = Jvx (x^yyy - ^{ g up y + g va y + g pa y}x^ 



+ ^pi^yy + g^yy + <ryT} 

o c \ 

' {g pv g pu + g pp g vu + g pa g vp } i, (4.36) 



AT 2 



j D \k p k v 2 k p 2 ki\ = Jvx (y p yy p y° + ^{<T<T + g w g va + g^g up } 



- ^p{g pu y p y a + gwyy + g^yy + g vp y»y° 

+ g ua y^y p + g pa yy}) i. (4.37) 



Note that these expressions are valid for arbitrary two-loop integrals and to use them all 
that needs to be done is to identify a, b, c, d p , e p and /, which can easily be read off from 
the graph, and to construct X p and y p . The powers of X{ and V can then be exchanged for 
scalar integrals with higher i/\ and higher D. This procedure is straightforward to implement 
in an algebraic program. 

Of course, rewriting the tensor integrals in this way is only useful if the relevant scalar 
integrals in higher D and with arbitrary powers of the propagators are known. For the 
pentabox, this is the case since repeated application of the integration-by-parts identi- 
ties ( |3T2|) and (|3.3|) can always be used to reduce the integrals to one of the master-topology 
integrals of Eqs. (|3.6| )- (|3~8| ), for which general results are given in the next section. 
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However, there is one further simplification of note. When in the numerator we have a 
dot product that can be rewritten in terms of the difference of propagators, it is easiest to 
directly cancel off one of the propagator powers. For example, if we have to deal with k\ - pi, 
it is more efficient to make the replacement 

2h- Pl = (h + Pl ) 2 -kl = A 2 -A 1 = 2~- I". (4.38) 

The integration-by-parts identities ( |3.2|) and ( |3.3|) can still be used to further reduce the 
powers of u%, 2/4, 2/5 and z/ 6 . 

Following this procedure, we can pinch completely the second propagator, so that we 
immediately obtain self-energy insertions to one-loop triangle graphs of the type shown in 
Fig. |2|. These can be written directly in terms of T functions 

J D (z/i,0, 2/3,0, 2/ 5 ,0, u 7 ; s,t) = n D (i/ 5 , 2/ 7 ) If ( 2/5 + v 7 - y,z/i,i/ 3 ;sV (4.39) 

2 

J D (0, 0,2/3,0, 2/5, u 6 ,u 7 ;s,t) = n D (2/ 3 ,2/ 7 ) If (2/5,2/3 + 2/7- y,f6;s) , (4.41) 
where the self-energy factor is given by 

u r U + V2 - f ) r (f - 2/i) r (f - ^) 

II D fa, 2/ 2 = -I t V , p \; n ^ V \ ^, (4.42) 

r (2/1) r (2/2) Y(D — v\ — 2/ 2 ) 

and the one-loop triangle integral with arbitrary powers and with one external mass scale 
Q 2 is given by (see for example |7|]) 

v y v y r (2/2) r (2/3) r (D - 2/i 23 ) 

where we have introduced the shorthand 2/^ = 2/j + 2/j, 2/^^ = i/j + 2/j + vj., etc. 



J D (0, 0,2/3, 2/ 4 ,0,2/ 6 , 2/7; s,t) = n D (2/ 6 , 2/7) Jg (2/4,2/6 + 2/7- ^,2/ 3 ;s) , (4.40) 



5 The master topologies 



In this section, we collect the explicit expressions in terms of hypergeometric functions for 
the general master-topology integrals A D (B D ) and C D , with arbitrary powers of propagators 
and general dimension D. We also show how to relate these to simpler pinched integrals 
and we compute explicitly the integrals that we need to evaluate all the tensor integrals in 
a closed analytical form and as an expansion in e (D = 4 — 2e). 
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Figure 2: One-loop self-energy insertion into a one-loop triangle. The propagators 
are labelled according to Eq. (2.S) and are each raised to the Vi power. 



5.1 Topology A 



The master-topology integral A D (z/i, z/ 2 , z/ 3 , is 5 , u 7 ; s, t) is shown in Fig. |] and, with respect 
to the propagator factors of Eq. ( [2.3|) , is defined as 

d D ki r d D k 2 



1 



in D/2 J m D/2 A" 1 A% 2 A3 3 A5 5 A" 7 ' 

This integral is related by a factor to the ordinary one-loop box integral 

A D {pi, v 2 , V3, vs, vi\ s, t) = U D (iy 5 , V 7 ) /f (u X , v 2 , Vz, v<j + v 7 - y ; s, t 



where IT (1/5, u 7 ) is given in Eq. ( [4.42 ) and 



if (ux, u 2 , u 3 , V±\ s,t) = J 



d u k x 



D/2 A u i A v i A u 3 A u i ' 
71 1 ^2 ^3 ^4 



ITi 



(5.i; 



(5.2) 



(5.3) 



Expressions for the one-loop box integral with general powers of the propagators are easily 
obtained ||, [?J |§. In the kinematic region \t\ < \s\ we find |§] 

r (^234 - f ) r (f - u 2M ) r (f - v X2A 



X 



+ 



r (z/ x ) r (z/ 3 ) rp- ^ 234 ) 

r, ( D 1 D 1 D f 

X 3 F 2 ( V2, "4, ^1234 - — , 1 + Z/234 ~ — , 1 + ^124 ~ ~, ~ 

r (z/234 - f ) r (f - z/ 34 ) r fa - ^) r (f - u 23 ) 

- Z/1234J 
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Figure 3: The master topology A. The propagators are labelled according to 



Eq. (2.S) and are each raised to the Uj power. 



D 



D 



D 



3^2 [Vi, — - y - ^23, 1 + ^1 - ^3, 1 + y - ^234, ~~ 



+ 



r (z/i24 - f ) r (f - u 14 ) r - v 3 ) r (f - v X2 
r (vj) r (u 2 ) r (z/ 4 ) r (d ~ ^1234 J 



X 



D 



3^2 ( ^3, y - ^14> y - ^12, 1 + V 3 -V 1 ,1 + — - Z/i 2 4, -- 



• (5-4) 



Expressions for If valid for \s\ < \t\ can be obtained by the exchanges V\ «-> z/ 2 , z/3 «-> 1/4 and 
s <-> £ in Eq. 



Although these integrals appear with a variety of integer values for U\, v 2 and z/3 in 
dimensions higher than D = 4 — 2e, we can systematically relate them to simpler integrals 
by the integration-by-parts identities. In fact, we find 



s v{l+lf = -(£>- z/i 2 334) if + + v 2 2 + + z/ 4 4+) 3~lf, 



(5.5) 



(D _ 2 - z/ 1223 ) z/ 2 2+/f = (D — 2 — u 13M ) v 4 4+lf + (y 2 - u 4 ) (u,l + + u 3 3+) if, (5.6) 
s u 3 3+lf = - (D - i/ 11234 ) if + (^2+ + z/ 3 3+ + */ 4 4+) l"if , (5.7) 



where we have used the shorthand notation If = If (fx, u 2 , v 3 , 2/4; s, t) and Vijjh = ^+2^'+^, 
etc. 

Repeated application of these identities reduces v\, v 2 and v 3 to unity together with 
integrals where one of the propagators is pinched out, yielding a one-loop triangle (see 

Eq. O) 

If (^i,^2,0, 1/4; = if (z/i,z/ 2 ,z/ 4 ;t) , (5.8) 
Jf (0, z/ 2 , 1/3, i/ 4 ; s, £) = if (1/3, z/ 2 , z/ 4 ; £) . (5.9) 
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Subsequent application of 



D 



(D - u 122U ) If + (^1+ + u 3 S+ + u 4 4+) 2~l 



- tD 

4 5 



can be used to control the power of z/ 4 and form the pinched triangle integral 

if (z/i, 0,z/ 3 ,^ 4 ; s,t) = (1/4,1^,1*3; s). 



(5.10) 



(5.11) 



Equation ( [5 . 1 0| ) should be used until z/ 4 = 2 — D/2, corresponding to 1/5 = Vj — 1, that is 
the master integral for the A topology. 



(5.12) 



The dimensional- shift identity for the master integral is obtained applying 

(z/ 3 3+ + z/ 2 2+ + v 1 l+)(vr 5+ + v 7 7+) + u 5 5+u 7 7 + 

to A D+2 (1, 1, 1, 1, 1; s, t), and reducing the right-hand side with the help of Eqs. ( ^.5|) -(^7| 
and ( PH : 



A D+2 {l,l,l,l,l;s,t) 



(£> - 4)st 2 



+ 



3(£> - 1)(3L> - 10)(3D - 8)(t + s 
s[{D-4)t+{2D-6)s) 



A D (1,1,1,1,1; s,t) 
D 



3(D -2)(D- 1)(3D -8)(t + s) 



+ 



D 



where 



3(D -A)(D -l)(t + s 



2--,l;f 1^(1,1) 



(5.13) 
(5.14) 



The master integral for the A topology has all propagators raised to unit powers and may 
written in terms of Gaussian hypergeometric functions as BJ] 



x 



I* (2 - f ) 


r 2 (f - 1) 


r(4- 


D)Y(D-3) 






r (d 


-2)r 







D-5 



r(f -2) 
r(2-f) 



.fl|i,i,£-i' + < 



+ (-s) D - 5 T (£> - 3) ^(1,2-^,^-1 



s + 1 



2 ' 2 s 

Making the expansion in e = -D/2 — 2 of the hypergeometric functions, we obtain 

r 3 (l-e)r(l + 2e) 



(5.15) 



A D (l,l,l,l,l;s,t) 



2 s (1 -2e) e 3 r(l - 3e) 
where A\ (s, t) and A 2 (s, t) are given respectively by: 



(s)- 2e A 1 (s,t) + (-t)- 2e A 2 (s,t) (5.16) 
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1) in the physical region s > 0, t < 0: 



A 1 (s,t) 



A 2 (s,t) 



1-e 2 



Li 2 ( — )-2( 2 



— e 



Li i ( ) +- S 2 2 



— e 

s + t 
t 



r . f 's + t\ _ fs + t 
Li 3 [— r- ) + Si, 2 



+ Si 



t / 
s + t 
t 



t 



+ 2C 2 Li 2 [' ^) - 9C 4 ] } + O (e 5 ) , 



s + t 



(5.17) 



1 + e log (-£) - e 2 Li 2 (^±*) - e 3 Li 3 (^) - e 4 Li, 
+ 0(e 5 ), (5.18) 



2) while in the region s < 0, t < 0: 



Ai (s,t) 



Si,: 



1 + e 2 
t 

s + t 



Lio 



S + t 



7T 

y 



— e 



2 Li, 



s + t 



Cs + log 



s + t 



Li 5 



+ -7T 

s + t/ 6 



+ e q 



Si, 



s + t 



2S 



2,2 



4v(i±i] + io g 



s + t 
s + t 



^ AT- ^ t 

+ 4Li 4 



23 



-7T 



t 



( 2L H s+t 



s + t/ 180 
t \ „ / t 



Si, : 



7T / t \ vr 2 2 /s + t\ I 

— Li 2 log 2 + - log 2 

3 Vs + t/ 4 B \ t J 2 & 



+ Y Ll 

+ o( e 5 r 



s + t 



Li- 



s + t 
t 



s + t 



(5.19) 



A 2 (s,t) = 1 + elog 



Lio 



s \ !, 

+ « log 

s + t/ 2 e 



s + t 



7T 

y 



— e 



+ e 4 



Li 3 



vr 2 /s + t 

H lot 

s + t 3 f 



1 , fs + t 
6 



LL 



7T 4 7T 2 2 /S + t 

log 2 

s + t 45 6 6 V s 



24 



log 4 



s + t 



+ e 



(5.20) 
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Figure 4: The master topology C. The propagators are labelled according to 



Eq. (2.S) and are each raised to the v\ power. 



5.2 Topology C 

The diagonal-box master-topology integral is shown in Fig. f| and with respect to the prop- 
agator factors of Eq. (|2.3|) is defined as 



d D h r d D k 2 



C D (V 2 , U 3 , "5, »6, "7] S,t) = f — f ~ 



l 



ztt d / 2 A? A? A? A? A?' 



(5.21) 



For this graph, it is easy to read off the values of Q and V of Eqs. (|4.8|) and (|4.10|) 

Q = x 3 x 6 x 7 s + x 2 x 5 x 7 1, (5.22) 

V = x 7 (x 2 + x 3 + x 5 + x 6 ) + (x 2 + x 3 )(x 5 + x 6 ). (5.23) 

The integral can be written as a single Mellin-Barnes integral (see for example Refs. |2|, |3j ) 



C D (u 2 ,u 3 ,u 5 ,u 6 ,u 7 ;s,t) 



r(z/ 2 )r(z/ 3 )r(i/ 5 )r(z/ 6 )r(z/ 7 ) 
r (f - u 7 ) r (f - z/ 56 ) r (f - v 23 ) ^ da 



T{D- v 237 ) //.-„;•,) I (ID - v 23567 ) J-'* -2k i 

xrp + i/j - z/23567 - a)T (z/ 235 67 - D + a) V (v 3 + a) T (v 6 + a) (^) . ( 5.2 I ) 



r (-«) r (D + V 2 - 1/23567 - a) 
t x a 



where the path of integration over a must be chosen so that to separate the poles coming 
from r (. . . — a) from those coming from T (. . . + a). In the kinematic region \t\ < \s\ the 
contour at the infinity must be closed to the right, and the expression we obtain is 

c» k „, „, , «) = (-!)» ^- r( f ^ 3 )r(f-, 56 )r( f -, ; ) 

r (z/ 7 ) r If D - 1/23567 
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+ 



r (^23567 -D)T(D- U 2367 ) T (D - Z/ 3567 ) 

r (z/ 2 ) r (z/ 5 ) r (d - z/ 567 ) r (d - z/ 237 ) 

X 3^2 (va, Va, ^23567 -D,l-D + U 23e7 , 1 - D + Z/ 3567 , -- 

r (^2 - vb) t(d- v 2%7 ) r (z/ 2367 - d) 



/f\D-U 236 7 / f 

X f - J 3^2 \ VS, D - U 267 , D - U 237 , 1 + D- Z/ 2367 , 1 + ffe ~ ^2, ~~ 

r (^ 5 - ^2) r (-D — z/ 357 ) r (^ 356 7 - d) 
r fa) r (i/ 5 ) r r p - */ 237 ) 

X f -J 3 -F 2 ( ^2, - Z/ 357 , D - Z/ 567 , 1 + D - Z/ 3567 , 1 + I/ 2 - Z/5, — 

(5.25) 

The solution valid when \s\ < \t\ can be obtained from Eq. (|5.25|) by the exchanges 

s t, z/ 2 <-> z/ 3 , z/ 5 <-> i/ 6 . (5.26) 

The expression for the diagonal box ( 5.25 ) has an apparent singularity when v 2 — i/ 5 is an 
integer which cancels in the actual evaluation of the diagram. 

In addition, although in principle we need to know this integral for arbitrary powers and 
arbitrary dimension, further simplifications occur. Integration by parts yields the following 
relations (C D = C D (z/ 2 , z/ 3 , u 5 , v$, v 7 \ s, t)) 

(D — 2 — 2v 23 ) u 2 2 + C D = (D — 2 — 2u 7 ) u 7 7 + C D - (D - 2 - 2v 23 ) v 3 3+C D , (5.27) 

(D-2- 2z/ 56 ) u 6 6 + C D = (D — 2 — 2v 7 ) u 7 7+C D - (D - 2 - 2z/ 56 ) z/ 5 5 + C D , (5.28) 

so that we can reduce both v 2 and z/ 6 to unity at the expense of increasing v 3 and z/ 5 together 
with v 7 . We now reduce u 3 and z/ 5 to unity using the relations 

s (D - 2 - 2v 23 ) u 3 3+C D = - (D - 1 - v 237 ) (3D - 2z/ 235667 ) C D (5.29) 

+ 2 (D - 1 - v 237 ) z/ 5 5+6-C D + (D — 2 — 2v 7 ) u 7 7 + 6~C D , 
t(D-2-2u m )u 5 5+C D = - (D - 1 - z/ 567 ) (3D - 2z/ 223567 ) C D (5.30) 

+ 2 (D - 1 - z/ 567 ) u 3 3 + 2~C D + (D — 2 — 2v 7 ) u 7 7 + 2~C D ', 
which, because i/ 2 and z/ 6 are already unity, produces simpler pinched integrals of the form 

C D (0, z/ 3 , u 5 , v 6 , v 7 \ s, t) = II D (z/ 3 , v 7 ) I 3 (v b , u 6 , v 3 + v 7 - y ; s 
C D (^2, ^3, ^5, 0, if; s, t) = Tl D (z/ 5 , z/ 7 ) (V 3 , IS 2 , v 5 + v 7 - ^; t ) , (5.31) 



where the self-energy factor II D is defined in Eq. (|4.42|) and the one-loop triangle I® is given 
by Eq. (TO - 



When the outer propagators have unit powers, we can reduce v 7 using 

st (D - 2 - 2v 7 ) u 7 7 + C D = -(s + t)(D-3-u 7 )(3D- 10 - 2z/ 7 ) C D 



(5.32) 



+ 2 (D - 3 - v 7 ) t5+6 - + s6+5 _ C 



iD 



+ (D — 2 — 2v 7 ) (tv 7 7 + 6- + sz/ 7 7+5 _ ) C 
This equation is only valid when v 2 — 1/3 — — i/ 6 = 1. 
Finally, applying 

d~ = z/ 7 7+ (z/ 2 2+ + u 3 3 + + u 5 5 + + z/ 6 6+) + (v 2 2+ + u 3 3+) (V 5 5 + + z/ 6 6 
to C D+2 (1, 1, 1, 1, 1; s, t), we can obtain the dimensional-shift equation 

(D - 4) W 



D 



(5.33) 



c D+2 (1,1,1,1,1; s,t) 



3(D - 3)(D - 2)(3D - 10)(3D - 8)(t + s) 



C"(l,l,l,l,l;s,t) 



+ 



s [(2D - 5)t + (D — 3)s] 
3(D-3)(D-2) 2 (t + S y 



D 



# 2-=-,l;a n»(l,l) 



t[(D-3)t + (2D-5)s} D y 
3(D-3)(D-2) 2 (t + s) 2 2 ' 



2 
D 



l;t \ 11" (1,1). (5.34) 



To summarize, we actually only need to know C D (1, 1, 1, 1, 1; s, t) rather than the more 
general case. With the propagator powers equal to unity, all of the 3 F 2 functions of Eq. ( |5.25| ) 
reduce to 2 Fi. To deal with the pole in (v 2 — z/ 5 ) we proceed as in Ref. J7j: we set v 2 = u 5 + 5, 
and, after performing an appropriate analytical continuation, we take the limit 5— >0. The 
final expression is given by 



C u (l,l,l,l,l;s,t) 



r(f -i)r(3-D)r 2 (f -2 



r(|D-5 



x 



t) D - 5 2 Fi [1,1, D -2, 



+ (s. 

If we make a series expansion in e, we obtain 

r 3 (l-e)r(l + 2e) 



< L '- 5 2^(1,1,^-2 



s + t 

t 

s + t 



(5.35) 



C u (l, 1,1, l,l;s,t) 



2 (s + t) T (1 - 3e) e 3 
where C (s,t) is given respectively by: 



(- a )-*C (*,*) + (-*)-*<?(*, a) 



(5.36) 
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1) in the physical region s > 0, t < 0: 



C(s,t) 



log ( -- I +2eLi 2 



s + t 



+ 4e 2 Li 3 ( —j- 



+ 8e 3 Li 4 



s + t 



+ O (e 4 ) 



2) while in the region s < 0, t < 0: 

CM) ' " 



log ( - ) - 2e 



+ 4e 2 



Li 3 



Li 4 

+ o( e + 



Li 2 

t 

s + t 
t 

s + t 



t 



s + t 

1 



6 



log d 



H lot 

24 6 



s + t 



t 

s + t 
t 



s + t 



7T 

y 



7T 



n 1 2 
U 



s + t 



t 

s + t 

t 



(5.37) 



7T 

45 



(5.38) 



Note that the prefactor of Eq. Q5.35 ) indicates that the integral diverges as 1/e 3 . However, 
the hypergeometric functions conspire to remove the leading divergence and we reproduce 
the result for C D (1, 1, 1, 1, 1; s, t) quoted in Ref. M. 



6 Conclusions 



Finally we summarize what we have accomplished in this paper. The pentabox with light- 
like external legs is one of the two-loop box graphs needed for the evaluation of Feynman 
diagrams for physical 2 — ^ 2 scattering processes, the others being the recently evaluated 
planar double-box graph |1|, 0] and non-planar double-box graph || as well as one-loop 
box integrals with bubble insertions on one of the propagators ||]. Using integration by 
parts identities, we have derived explicit analytic expressions for the scalar pentabox graph 
J D (1, 1, 1, 1, 1, 1, 1; s, t) in terms of two simpler box integrals (A and C) with fewer propaga- 
tors and given an explicit expansion in e in terms of polylogarithm functions (see Eq. fl3.13| )). 

In Section [| we demonstrated how generic two-loop integrals could be expressed directly 
in terms of integrals in higher dimension with shifted propagator powers. Similar results 
have been given by Tarasov |J] using differential operators; however we believe our method 
based on completing the square and Gaussian integration to be more direct and simpler 
to apply to generic integrals where the quantities a, b, c, d^, and / (see Eq. (|4.3|) ) can 
easily be read off from the graph and used to construct the vectors X^ 1 and y u (defined in 
Eq. (|4.7|) ), which therefore are linear in both the external momenta and the As detailed 
in Section |], the powers of Xi and V present in Eqs. (4.22) and (4.25)-(4.37) can then be 
straightforwardly exchanged for scalar integrals with higher V{ and higher D. 
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In the case of the pentabox, the tensor integrals can also be written in terms of the same 
master-topology integrals as the scalar pentabox, but with different powers of propagators 
and in higher dimension. We give quite general analytic formulae for the master-topology 
integrals with arbitrary dimension and propagators. For the A topology we quote the results 
of 0] while for the C topology we employ the Mellin-Barnes integral representation. Other 
methods would also suffice. In fact, it is not necessary to know the master integrals in 
their full generality since further integration-by-parts identities relate the various integrals. 
Explicit expressions for topology A and C are given in Eqs. ( |5.15|) and (|5.35|) respectively, 



together with instructions on how to apply them to more general integrals. 

Finally, we note that the recurrence relations necessary for the complete reduction of 
tensor integrals for the planar-box graph were described in Ref. 0. However, the situation 
for the non-planar box is less well developed. At present only one master integral is known, 
and the solution of the recurrence relations for that topology is an open problem. 
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A Poly logarithms 

The purpose of this appendix is to define the generalised polylogarithms that occur in the 
expansion in e of the pentabox scalar and tensor loop integrals and to give useful identities 



amongst the polylogarithms. In Section |A.1| we give the definitions of the polylogarithm 
functions S ntP (x). These functions are real when x < 1 but they develop an imaginary part 
for x > 1. Analytic continuation formulae are given in Section A.2 . Finally, useful identities 
between polylogarithms are listed in Section [A. 3. 



A.l Definition 



The generalised polylogarithms of Nielsen are defined by 



Sn, P {x =7 -T7-T / at , n,p>l, x<l. A.l 

(n — ly.p] Jo t 
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For p = 1 we find the usual polylogarithms 

Sn-1,1 (a;) = Li n (x). 

The Snap's with argument x, 1 — x and 1/x can be related to each other via |§ 



S n , P (1 x ) — E ~7j 



+ 



s=0 

i-iy 

n! p! 



p-1 /jNr 

S n _ SiP (l)-£^-^log r (a;) S p _ 



r=0 



log n (l-x)log p (a;) 



n i p 



:-iru-iy± { -^io g r (-x)( n+s r 1 ) s n+ ,_, 



s=0 r=0 



s — r 



+ E L -^— ^(-x) C n _ PiP + 7^rrlog n+p (-n 



r=0 



T ! 



(n + p)\ 



with 



p-i 



n, p 



(-ir +1 E(-ir r ( n+ r r 1 ) s„ + , P _ r (-i) 

+ -(_!)») S n , p (-1). 

The usual binomial coefficient is defined as 



V 6 / 6! (a -6)! 
We also need the definition of the Riemann Zeta functions 

OO 1 



8=1 



and in particular 



O 



C 3 = 1.20206... C4 = — 

90 



A. 2 Analytic continuation formulae 



For x > 1, the following analytic continuations should be used 
Li 2 (x + i0) = -Li 2 Q-Ilog 2 (a;) + y + mlog(a;) 

/ 1 \ 1 7T 2 Z7T 

Li 3 (ar + iO) = Li 3 ^-J - - log 3 (x) + y log(x) + y log 2 (x) 
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Li 4 (x + iO) 
Si, 2 (x + i0) 



S 2 ,2(x + iO) 



u 1 'J) " ki log4(;r) + y log2(;r) + h + j log3{x) 



(A.ll) 



Si. o ( -) + Li 3 + log(x)Li 2 Q 



1 7T 2 

glog 3 (x)- — l g(x) + C 3 



+ i7T 



7T 



— - Li 2 ( - 

o \x 



1\ 1 



-log (x) 



(A.12) 



Si, 3 (x + <0) = -S^aQ+S^Q+log^S^Q-L^Q-log^LiaQ 

1, 2/ nt- f l \ T 2 T - 1 , 4 / \ 7T 2 , 2/ \ 197T 4 



+ Z7T 



(A.13) 



7T 



+ Cs l0g(x) + — + Z7T 



Lb 



x 



1 7T 2 

-log 3 (x) + y log(x) - (3 



(A.14) 



A. 3 Useful identities 



For the expansions of hypergeometric functions given in this paper, the argument of the 
polylogarithms is always either (s + t)/s or (s + t)/t, which are related by the following 
identities 



Li 2 
Li 3 
Li 4 



x 



Si, 2 
Si, 3 
S2,2 



X- 1 

X 

x — 1 

X 

x — 1 



X 

x-1 

X 

x- 1 

' X 
,x — 1 



- Li 2 (x) - ^ log 2 (1 - x) , 



1 



-Li 3 (x) + Si, 2 (x)+log(l-x) Li 2 (x) + ^ log 3 (1-x), 

o 



(A.15) 
(A.16) 



Li 4 (x) + S 2 , 2 (x) - Si, 3 (x) + log (1-x) Li 3 (x) -log(l - x) Si, 2 (x) 

(A.17) 



-\ log 2 (1-x) Li 2 (x)-llog 4 (l-x), 



Si, 2 (x)-ilog 3 (l-x), 
-Si, 3 W-^log 4 (l-x), 



1 



S 2 , 2 (x) - 2 Si, 3 (x) - log (1 - x) Si, 2 (x) + ^ log 4 (1 - x) 



(A.18) 
(A.19) 
(A.20) 
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